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In the limit of infinite yield time for stresses, the hydrody- 
namic equations for viscoelastic, Non-Newtonian liquids such 
as polymer melts must reduce to that for solids. This piece 
of information suffices to uniquely determine the nonlinear 
convective derivative, an ongoing point of contention in the 
rheology literature. 
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Hydrodynamics, the macroscopic description of con- 
densed systems in the low frequency, long wavelength 
limit |l|,|[, is a well-established, systematic approach in 
condensed matter physics. This method has been applied 
over the past decades to many systems, including simple 
fluids ||, various liquid crystals pp|, and the superfluid 
phases of 3 He and 4 He ||. 

All viscoelastic non-Newtonian fluids behave as New- 
tonian ones at low frequencies, and as solids at higher fre- 
quencies. A consistent hydrodynamic description needs 
to reflect this fact and must therefore contain, as spe- 
cial cases, both the hydrodynamic theory for isotropic 
liquids and solids. The liquid limit is well heeded in the 
polymer literature and universally correctly implemented 
HQ. The solid limit is problematic, as we shall see, and 
compatibility especially in the nonlinear regime of large 
displacements and rotations has so far proven elusive. 
The reason behind it is probably the lack of a consistent 
hydrodynamic theory for solids. 

The last sentence may come as a surprise, but the point 
we are making here is: Although both the nonlinear elas- 
ticity theory || and the linearized hydrodynamics for 
crystals |l],|8|Q are well known and established, a consis- 
tent hydrodynamic theory that includes both nonlinear 
and irreversible terms is not - in spite of some insightful 
papers |l(|. One of the obstacles is that such a the- 
ory necessarily employs a strain tensor different from the 
one customarily used The usual strain tensor is of 
the Lagrange type, derived from equations of motion for 
mass points, while a framework to set up hydrodynamic 
equations including dissipative terms only exists in the 
Eulerian description - which considers evolution of field 
variables at spatial points. Consistency forbids a mix- 
ing of both descriptions and requires an Eulerian strain 
tensor jllj . (We note that the linear hydrodynamic the- 
ory may mix both descriptions, as the smallness of the 



displacements ensures that the discrepancy is negligible.) 

The presentation of the nonlinear hydrodynamic the- 
ory for solids is what we shall do first. Then these equa- 
tions are generalized for non-Newtonian fluids by adding 
relaxation-type terms to account for a finite yield time 
of the stresses, such that in the high frequency limit the 
theory is unchanged, but in the low frequency limit only 
the terms of the isotropic liquid hydrodynamics remain. 
So, by ensuring the valid liquid and solid limits, this 
approach leads to the correct hydrodynamic theory for 
any liquids displaying viscoelasticity. ft is of great im- 
portance for rheology, as many competing theories exist, 
which differ especially in their respective nonlinear con- 
vective derivatives. All are employed in the literature, 
with scant hope for consensus. 

Comparing our result with the literature, we find that 
none of the convective nonlinearities suggested in the 
rheology literature ||f7| contains the correct solid- limit, 
putting justified doubts on their validity. (We do obtain, 
for the Eulerian strain and in the limit in which it is 
small, the so called "upper convected derivative".) 

Note that the insistence on the solid limit also deter- 
mines the choice of variables, being that of a solid: the 
strain and the conserved quantities. As the solid is the 
more complicated of the two limits, there is no reason 
to, in addition, take the stress as an independent vari- 
able, as most theories in the rheology literature do, see 
for instance Chap. 7-9 in the first of Aside from un- 
necessarily making a derived quantity independent, this 
approach also commits a cardinal hydrodynamic sin, be- 
cause the stress as a hydrodynamic flux contains reactive 
and dissipative parts, and does not possess a well-defined 
time reversal parity - without which we have no way of 
applying the Onsager relations. 

We now introduce the nonlinear hydrodynamic theory 
of solids. A proper description relies on two coordinates: 
the actual spatial coordinate n, specifying a point in an 
elastic body, and the coordinate this point possesses in 
the absence of any stresses. More carefully, starting from 
a stress-free elastic body, we consider a point with the 
initial coordinate a^. As the body is translated, rotated, 
compressed and sheared, this point is displaced to r% - 
especially in soft matter generally rather remote from 
ai. Since all points of the body have a unique pair of a% 
and rj, the function ri(a m ) is unique and invertible, the 
result of which is denoted as <2j(r m ). For briefness, we 
shall refer to all rj as the real space, and to all <Zj as the 
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initial space. 

As discussed in most books on elasticity theory, see 
eg ||, the elastic energy depends on the change in the 
distance between any two neighbouring points, from daf 
to drf. Defining the displacement vector as Ui(a m ) — 
fi{a m ) — cii, and the strain tensor as Ufa — \\duijdak + 
diik/dcii + (duj/dai) ■ (duj/ddk)], we have dr?(a m ) — 
daf = 2[/^,daidafe, and know to take the energy den- 
sity as a function of U^ k , to lowest order simply as 
e = ^KikjmU^Uj^. The important point here is that 
both the strain tensor and the energy density are func- 
tions of the initial coordinate a m , a notation that we shall 
refer to as Lagrangian - hence the superscript in U^ k . 

Contrast this with the energy density of an isotropic 
liquid in its rest frame, a function of the mass and entropy 
density, e(p, s) - or equivalently, de — Tds + pdp. All 
variables, including the conjugate ones, temperature T 
and chemical potential fi, are here functions of the real 
coordinate r m . As a result, the spatial dependence of 
(say) the temperature is quite independent of the liquid's 
compressional state. This is the Euler notation, and its 
basic advantage is that physics, which we insist must 
be local, is also expressed in local terms, accounted for 
by quantities at the real coordinates r m . Consider for 
instance the diffusive heat current, which is given by the 
local gradient of the temperature, ~ dT(r m ) / dr k , only 
in the Eulerian description. 

Returning to solids, we have two choices: First, take 
all variables including especially the temperature and 
chemical potential as functions of a m , and employ them 
with the strain tensor l/R. This would be consistent, 
but highly inconvenient. For instance, the heat cur- 
rent ~ dT(r m ) I 'dri at the real space point r m now 
presumes the knowledge (not usually available) of the 
global transformation, r m <-» a m , as dT(r m )/dn = 
[dT(a m )/dak](dak/dri). Similarly, with g the momen- 
tum density, the angular momentum density is r(a m ) x 
g(a m ) rather than a x g(a m ). (If the system is only 
weakly deformed, with it, = — a$ small, the above dif- 
ferences between and ai may be neglected to linear 
order.) Finally, more specific to the issue at hand, our 
equations need to contain both the elasticity theory and 
the liquid hydrodynamics. The latter, however, is usually 
and concisely given in the Euler notation. 

The second, and the only actually viable, choice is 
to take all variables including the strain tensor in the 
local, Eulerian notation, as functions of r m . We shall 
therefore employ the Eulerian strain tensor |lljl2|, in- 
troduced via drf — daf(r m ) = 2Uik(r m )dridrk, where 
U ik (r m ) = \[dui/drk + du k /dn - {du J /dr i ){du j /dr k )] 
and Ui(r m ) = n - Oi(r m ). 

In this context, there is a second, somewhat subtle 
point: We need to eliminate the displacement field Ui, 
and deal exclusively with the initial coordinate ai(r) 
and the strain Uik in the elasticity theory. This is 
possible because starting again from drf — daf(r m ) = 
2Uik(r m )dridrk, we find that the strain tensor may be 



written as 



U ik = ^[5ik - (da a /dr k )(da a /dri)}, 



(1) 



with no need whatever for a detour via Ui. This is neces- 
sary because the introduction of Ui destroys a symmetry 
and represents an arbitrary (gauge) choice. As discussed, 
di and n are vectors of different spaces, so they trans- 
form as vectors under rotations in initial and real space, 
respectively. The introduction of the displacement fixes 
both spaces with respect to each other, and prohibits the 
rotation of either space alone. On the other hand, the 
elastic energy is independent of the orientation of the 
initial space, the fictitious unstressed body. Given any 
transformation a m «-> r m , we should still be free to take 
a global but arbitrary rotation of all a^, ie, rotate the 
initial space with respect to the real space. Therefore, 
a a and rj are indeed vectors of two different spaces, and 
a quantity such as Vfca a = da a /dri is a vector both in 
real and initial space, (a bi-vector,) and not a second 
rank tensor. We use Latin and Greek indices to denote 
the components (x,y,z) in real space, and (1,2,3) in initial 
space, respectively. (Clearly, this renders the fact that 
the displacement — a a is an oxymoron rather obvious.) 

We now proceed to understand that the bi-vector V \a a 
not only contains the information about the strain, as 
shown in Eq (|l|), but also that about the local orientation. 
The polar decomposition theorem (cf W. Noll, p. 65 ff, 
Vol.2 of §) states 



— Raj ^ij 1 



(2) 



where R a j is the rotation matrix that rotates the local 
preferred directions in real space back to the global ones 
in initial space, while is a symmetric matrix that de- 
viates from Sij only for finite strains. Consider first the 



unstrained case 5. 



Because of da a — (Via a )dr.i 



with da/^ — drf, the gradient Via Q is indeed a ro- 
tation matrix R a j, and must satisfy R a jR a k = Sjk, 



Raj Rfij — 

of which is 



6 a /3. For finite strains, Eq. ([T]) implies 
R ak E ik R ai Eij = HifcSfej, the square root 



Eij = 5 i: j - 2Uij w (Sij - Ui-j - \U lk U k j ■ ■ •)• (3) 

[This expansion is valid for small strains C/y, but ar- 
bitrary rotations R a j- The square root of a matrix is 
defined by its power series. One can verify Eq (^) by 
calculating S ifc S fcj -.] 

In accounting for solid behaviour, we need to keep 
track of the local preferred directions, or R a i, which may 
vary considerably by accumulation over a long distance, 
even if the strain is small - think of a sheet of single 
crystal, slightly bent over a long stretch to form a tube 
of large radius. Let us consider as an example the har- 
monic approximation for the energy E = J e dV, 

£ = —KijkmUijUkm = \Kaf3~jsUapU~f8 (4) 
Uij = RaiR/3jU a p, Kikjra = RaiRpjRykRSmKapyS (5) 
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where U a p and K a p-fS are the attendant quantities in 
the initial space, (e has the same form in both spaces be- 
cause R a i annihilates pairwise.) The elements of K a p 7 s 
are constant. A cubic crystal for instance has three in- 
dependent elements, of which one is given as Kim — 
K2222 = -K3333, implying that the compressional energy 
is the same along the three initial space directions 1, 2, 
and 3. This is not the case for a hexagonal crystal, for 
which the compressional energy along 3 is different from 
1 or 2. The real space matrix Kij km depends on R a i and 
varies in space, as Eq(||) shows, because the symmetry 
axis 3 (of a hexagonal crystal) may in real space point 
in any direction, and vary spatially. So e is a function of 
Uij and R a i. Writing de = ^ydZ/y + XaidR a i, the conju- 
gate variables ^y = Ki jkm Ukm and Xai = 2RpjFL k Rs m 
Ka/s-ygUijUkm are given by differentiating Eqs (|J) and 
(||). The energy e depends on Uij and R a % beyond the 
validity of Eq (||), so de = ijdllij + \aidRai is generally 
valid - though the explicit form of and Xai will vary. 

The term XaidRai is rarely included in the energy 
for solids ||, which renders the resultant formulas valid 
only for small deformations or isotropic solids. For the 
latter systems there is no preferred direction to keep 
track of locally, so Kijkm will not depend on R a t, and 
we can set Xai = 0. This is easiest seen in the har- 
monic approximation, Eq (Q), where K a f3jS — (Kl — 
K T /3)5 a/ 35 1 s + K T /2(S a ^Sf3s + SasSp^) due to isotropy. 
Inserting this into Eq (f5|), we again obtain Kijkm — 
(K L - K T /?,)8 l j5km + K T /2(8ikS jm + S im Sjk), manifestly 
independent of R a i. 

Returning to anisotropic systems, the 9 variables of 
Vid Q are equivalent to the 3 of R a i and the 6 of 
Uij (or Sjj), see Eqs. (||) and (J3J) . So we can con- 
veniently write ^JijdUij + XaidRai = ^aidVjCta where 
ipai = ^km{dUkm/dVia a ) + xpk{dRpk/dW l a a ). Under 
a real space rotation of the angle d9i, scalars are in- 
variant, de = ijj a i&Via a = 0, but vectors and tensors 
are not, dVja a = e ijk V a a d6 k , so ipai^ijk^ j® a — 0, or 
(ipajVida) = (i <-> j). Similarly, dR m = e ijk RajdO kl 
dU m = eijkUj m d9k + e m jkUijd9k, so (^i k Uj k + ^ k iU k j + 
XaiRaj) = (i <-> j)- (These constraints on ip a j, ^?U, and 
Xai can be used to show the symmetry of the stress tensor 
<Jij below.) 

More generally, e also depends on the mass, entropy 
and momentum density, p, s, and So the final ther- 
modynamic expression for an elastic medium is 

de = Tds + pdp + Vidgi + ^crfdV \a a - (6) 

Turning now to dynamics, the equation of motion for 
a a is 

^a a =a a + v k V k a a = ~Y a . (7) 

In equilibrium, with the dissipative contribution Y a van- 
ishing, this equation simply states the fact that the initial 
coordinate a a of a mass point does not change when one 
moves with it. 



The entropy production s + Vi(svi — fa) — R/T, con- 
servation of mass and momentum, p + Vi(pvi) = 0, 
cji + Vj(er.y ; — afj) = 0, and Eq. (fjj) represent the complete 
hydrodynamic theory of solids, where 

Uij = pSij + Vig-j + IpajVida, (8) 

R = /<V<T + afjAij - YaV k ^ak, (9) 

[with A ik = |(Vi«fe+VfeUi), P = -e + Ts + p,p + Vigi] are 
unambiguously given by thermodynamics, Eq. (o), via 
the hydrodynamic standard procedure. Eq. (^|) implies 
that the three fluxes fi,Oy,Y a are linear combinations 
of the three forces VjT, Ay, Vfc-0afc- These give rise, re- 
spectively, to the dissipative phenomena of diffusive heat 
current, viscous stress, and defect diffusion. The struc- 
ture of the linear combination, ie the symmetry of the 
Onsager coefficients, are given by the symmetry group of 
the crystal ||. For isotropic solids, we have only diag- 
onal terms, especially /j ~ V^T and Y a ~ V k ip a k- As 
discussed at length in , it is incorrect to take the lat- 
ter contribution as zero: The initial coordinate of a mass 
point may change, a a =/= 0, in the absence of any mass 
current, Vi — 0, when there is diffusive motion of va- 
cancies. Conversely, motions of interstitials involve mass 
current, «j ^ 0, but no change of crystal points, a a = 0. 

Since Via a as a variable is completely equivalent to Uij 
and Rai 7 the equation of motion (0) for d a may always 
be rewritten as two equations of motion, for Uij and R a i- 
Though rather more complicated, this is certainly closer 
to the conventional elasticity theory. With the help of 
Eqs. (0) and (0), we rewrite Eq. (0) as 

2f t U i0 = [EjiEikViVh + R ak E jk ViY a ] + [*<-» j], (10) 

Ral^Raj — ^Ij^i^j i? a ;V \Y a dt L ~ J ^' (H) 

which may be approximated by taking Uij , Aij , Y a and 
Raj 37 R a i — ^ij as small quantities (with 2wy = VjVi — 
ViVj). To second order in the small quantities, though 
neglecting terms of order U k jViY a , the result is 

&Uij - A tj - [±(ViY a )R aj - [ViV k )U kj ] + (12) 

Raj^Rai - Uij = [\R a i^ 3 Ya + \U j k A ik \ - [i<-*j] (13) 

Written in the conjugate variables of Uij and R a i, the 
stress tensor Eq. (||) reads 

cry = pSij + !■,</, - + ,/ + ^kjUik 

+XajRai + \Xak{U ki R a j + U k] Rai). (14) 

This ends the presentation of the hydrodynamic the- 
ory of solids. The noteworthy point is: The derivation 
is completely cogent, as not a single step in it is dis- 
cretionary; hence the above set of differential equations, 
given between Eq (Q) and (||) , is unique - any other the- 
ory is either algebraically equivalent, or wrong. Con- 
versely, these equations account for any solid system, in- 
cluding crystals of all symmetry groups and glasses. This 
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pertains especially to the nonlinear structure, important 
if one is to account for large displacement and rotation, 
strong compression and shear. These are usually small 
in bulk crystals, but quite large in complex liquids. In 
awareness of this, many nonlinear models for convective- 
like nonlinearities have been suggested , though none 
was constructed to contain the nonlinear solid limit. 

To generalize our results to visco-elastic Non- 
Newtonian fluids, we note first that the solid hydrody- 
namics contains the liquid hydrodynamics, and one can 
reduce the former to the latter by setting to zero the 
elastic stress ip a i — de/d(Via a ). Confining ourselves 
to isotropic systems, it suffices to set $y = 0, because 
Xai = already holds. Taking Uij as the variable that 
relaxes as long as ^ij ^ 0, we connect the isotropic solid 
dynamics to fluid dynamics such that the former holds in 
the high frequency regime (where the relaxation is negli- 
gible) and the latter in the low frequency regime (where 
relaxation is dominant). Therefore, we proceed by allow- 
ing a relaxation term Xij in Eq. (|l^) , Uij + • ■ ■ = X^ . 
It leads to an additional term in the entropy production, 

Eq (§), R = A, .>!',,. which implies X = -£ a^ 1 

in an expansion, or to lowest order, = -ax^Jj and 
Xkk = —oil® u- (The superscript denotes the traceless 
part of the given tensor.) So we have 

&Uij - Aij + [(ViV k )U k j - ±(V t Y a )R aj 

= -a T ^ij - ai*fefc%/3, (15) 

with aT,aL denoting two transport coefficients. To un- 
derstand the added terms, one can use the example of 
the harmonic approximation, Eq (Q), yielding ax^ij = 
cx T K T U% = U^/tt and a L V kk = a L K L U kk = U kk /r L . 
Clearly, this implies relaxation for L7y and U kk , with 
the respective relaxation times tt and T£. (In prin- 
ciple, there are two thermodynamic cross derivatives, 
S«£ kk = K p 5p + K T 5T.) 

Note the universality of the results, especially the con- 
vective terms ~ (VjUfc), which remarkably are not pre- 
ceded by any material-dependent coefficients. Their form 
is independent from the above approximation for X^ and 
will remain unchanged even if additional variables are in- 
troduced, eg to account for the material-dependent rheo- 
logical behavior such as shear thinning and normal stress 
differences. 

To account for large deformation, rotation and velocity, 
many different nonlinearities, as mentioned, have been 
suggested and implemented in the rheology literature, of 
which the two more popular ones are the upper and lower 
convective derivatives. Denoting an arbitrary matrix as 
(*), the former is defined as D u (*) = (d/dt+v • V)(*) + 
(Vw)(*) + (*)(Vv) T , the latter as D e (*) = (d/dt + 
v • V)(*) — (Vi>) T (*) — (*)(Vt>). Both are derived by 
invoking some variant of a postulated general principle, 
usually referred to as the "material frame independence" . 
In the rheology literature (*) is the stress tensor, 

taken as independent, but in principle it could also be 
the strain tensor. 



Reviewing the many equations of motion considered 
above, it is easy to see that Eqs. (|l^) and (|l^) for 
can indeed be written as D u U — A = 0(VY, \I/). None 
of the other equations may be brought into this form: It 
is not valid for a a and R a i, see Eqs. (fj], [TT], |lj|); nor for 
the exact equations Uij , Eq. ( |Io|) . Especially, it does not 
hold for the stress tensor. 

This is a rather serious shortcoming and subjects all 
those descriptions that include upper and lower con- 
verted derivatives, combinations thereof, or other kinds 
of quadratic nonlinearities to grave doubts. To overcome 
these, the authors really need to convincingly argue why 
their postulated general principle overrules the simple 
and physical requirement that, for infinite yield time of 
the stress, the dynamics of non-Newtonian liquids such 
as polymer melts is that of an isotropic elastic medium. 
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